Critical exponents and equation of state of three-dimensional spin models 
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Three-dimensional spin models of the Ising and XY universality classes are studied by a combination of high- 
temperature expansions and Monte Carlo simulations. Critical exponents are determined to very high precision. 
Scaling amplitude ratios are computed via the critical equation of state. Our results are compared with other 
theoretical computations and with experiments, with special emphasis on the A transition of 4 He. 



1. INTRODUCTION 

The notion of universality is central to the 
modern understanding of critical phenomena. It 
is therefore very important to compare high- 
precision theoretical and experimental determi- 
nations of universal quantities, such as critical 
exponents or universal ratios of amplitudes, for 
systems belonging to the same universality class. 

Critical exponents and amplitudes parametrize 
the singular behavior of thermodynamical quan- 
tities in the vicinity of a critical point. In 
the high-temperature (symmetric) phase t > 0, 
C H « A+\t\~ a , X « C+\t\-\ $ a f+\t\- v , where 
t = (T — T c )/T c is the reduced temperature, Ch 
is the specific heat, \ is the magnetic suscepti- 
bility, and £ is the correlation length. In the 
low-temperature (broken) phase t < 0, H — ► 0, 
C H « A-\t\- a , X « C-\t\-i, £ « f-\t\-" (in the 
case of Ising), M w -B|t|^, where M is the magne- 
tization. On the critical isotherm t = 0, H ^ 0, 
X « C c \H\-~<l f3S , £ « r\H\- v l< 3& . At the critical 
point t = 0, H = at nonzero momentum, the 
two-point function behaves like G(q) « Dq v ~ 2 . 

The critical exponents are universal, and are in- 
dependent of the phase; they are related by the 
scaling and hyperscaling relations. The ampli- 
tudes are not universal, and their value depends 
on the phase; it is however possible to define uni- 
versal ratios of amplitudes, which are indepen- 
dent of the normalization of H, M, and T. 

The universality classes of Ising and XY in 



three dimensions have been the subject of many 
theoretical studies. Nonetheless, we believe that 
further refinement is worthwile: many critical 
phenomena in nature fall into these classes, and 
the precision of experiments is ever improving; 
moreover, several theoretical techniques can be 
applied and compared to each other. 

High-temperature (HT) series expansion is one 
of the oldest and most successful approaches to 
the study of critical phenomena. We are extend- 
ing the length of the series available for wide fam- 
ilies of models belonging to the classes of univer- 
sality we are interested in; so far we computed the 
two-point function of the three-dimensional Ising 
class to 25th order on the bec lattice, and four-, 
six-... point functions to 21st, 19th... order jt). 
Work is in progress on the sc lattice and on the 
XY class. 

The precision of the results which can be ex- 
tracted from long HT series is mainly limited by 
the presence of confluent corrections with nonin- 
teger exponents. Let us consider, e.g., the mag- 
netic susceptibility %; near the critical tempera- 
ture, it behaves like 

X = Ct-~<(1 + a ,it + ao^t 2 + ... 

+ ai,i* A + ai, 2 t 2A + ... + a 2A t A2 + ...). 

While the exponents A, A 2 , ... are universal, the 
coefficients a are model dependent. For the mod- 
els we are interested in, A ~ 0.5 and A 2 > 2A; 
therefore it is very helpful to select one-parameter 
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families of models, and tune the irrelevant param- 
eter A to the special value A* for which ai.i = 0; 
we will call such models "improved" . 

Monte Carlo (MC) algorithms and finite-size 
scaling techniques are very effective in the deter- 
mination of A* and (3 C , but not as effective in the 
computation of critical exponents or other univer- 
sal quantities. On the other hand, the analysis of 
HT series is very effective in computing universal 
quantities, but not in computing A* and (3 C . 

The strength of the two methods can be com- 
bined by computing A* and j3 c by MC, and feed- 
ing the resulting values into the analysis of HT 
series (by "biasing" the analysis); this greatly im- 
proves the quality of the results. 

In order to keep systematic errors under con- 
trol, we always select several different families of 
models in the same universality class and check 
that they give compatible results for universal 
quantities. 

2. CRITICAL EXPONENTS 

Without further discussion, we present in Ta- 
ble [j] a selection of results for the critical expo- 
nents 7, v, and r\ of the three-dimensional Ising 
model; for other exponents, see Rcf. §. We com- 
pare the most precise theoretical results and ex- 
periments. 1HT denotes our results |Q; HT is 
a "traditional" HT determination ||; MC are 
Monte Carlo results Q; FT are results from a 
g expansion in fixed dimension |5j . Experimental 
results are LV for liquid-vapor transitions; BM for 
binary mixtures; MS for uniaxial magnetic sys- 
tems; MI for micellar systems; cf. Rcfs. || and 
Q for bibliographical details. The agreement be- 
tween the different determinations is overall sat- 
isfactory. 

On the theoretical side, similar techniques can 
be applied to the XY model, with results of com- 
parable quality. The experimental situation is 
quite different: one extremely precise experiment 
on the A transition of 4 He @ overshadows the 
field. We present results for the critical expo- 
nents 7, ?7, and a (we remind that dv = 2 — a) in 
Table || (cf. footnote 2 in Ref. ^ for discussion of 
the experimental results). Theoretical results are 
taken from Refs. § (IHT), [§ (IHT*), |] (HT), 



Table 1 

Comparison of determinations of critical expo- 
nents of the three-dimensional Ising model. 





7 


V 


V 


IHT 


1.2371(4) 


0.63002(23) 


0.0364(4) 


HT 


1.2384(6) 


0.6308(5) 




MC 


1.2367(11) 


0.6296(7) 


0.0358(9) 


FT 


1.2405(15) 


0.6300(15) 


0.032(3) 


LV 


1.233(10) 




0.042(6) 


BM 


1.228(39) 


0.628(8) 


0.0300(15) 


MS 


1.25(2) 


0.64(1) 




MI 


1.237(7) 


0.630(12) 


0.039(4) 


Table 2 








Comparison of determinations of critical expo- 


nents of the three-dimensional XY model. 




7 


V 


ot 


4 Hc 






— 0.01056(38) 


IHT 


1.3179(11) 


0.0381(3) 


— 0.0150(17) 


IHT* 


1.3177(5) 


0.0380(4) 


-0.0146(8) 


HT 


1.322(3) 


0.039(7) 


-0.022(6) 


MC 


1.3177(10) 


0.0380(5) 


-0.0148(15) 


FT 


1.3169(20) 


0.0354(25) 


-0.011(4) 



H (MC), and f§ (FT). There is disagreement be- 
tween IHT* and experiment; it would be interest- 
ing to improve further the theoretical computa- 
tion, and to have an independent confirmation of 
the experimental measurement. 

3. CRITICAL EQUATION OF STATE 

The critical equation of state is a relation be- 
tween thermodynamical quantities which is valid 
in both phases in the neighborhood of the critical 
point (cf., e.g., Ref. 0). 

In order to determine the critical equation 
of state, we start from the effective potential 
(Helmholtz free energy) 

T(M) = ME - ^ log Z{H). 

In the high-temperature phase, J- can be ex- 
panded in powers of M 2 around M = 0. By 
choosing appropriately the normalizations of the 
renormalized quantities, and using the "second 
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moment" mass m as mass scale, we can write 
Af = f(M)-F(0) = — A(z), 

94, 



3 >3 K J > 



,2/ 



where z is the (rescaled) zero-momentum vac- 
uum expectation value of the renormalized field, 
r 2j = 92j/gi~ 1 , and g^j is the renormalized 
zero-momentum 2j-point coupling constant. The 
(universal) critical limit of 54 and r 2 j can be 
computed from the HT expansion of the zero- 
momentum 2j-point Green's function; for the 
Ising model, we obtain [|| 

.94 = 23.49(4), r 6 = 2.048(5), 
r 8 = 2.28(8), no = -13(4). 

The equation of state can now be written as 



zocMr' 3 ; 



(1) 



The analyticity properties of F(z) are constrained 
by Griffiths' analyticity conditions on H(M,t). 

It is possible to implement all analyticity and 
scaling properties of the critical equation of state 



introducing a parametric representation 1 1 



M = m o R 9, 

t = R(l- 6> 2 ), 
H = h R f3S h(9), 



h(9) 



0{9 3 



The following correspondences should be noticed: 

= — > t >0,M = 0; 
9 = 1 — > t = 0; 
9 = 9 Q — > t < 0,M = M , 

where 9q is the first positive zero of h{9). The 
analytic properties of the equation of state are re- 
produced if h(9) is analytic in the interval [0, 9o\. 

Combining the parametric representation with 
Eq. (Q), we obtain 

z = P 9 2 {\-9 2 )- 13 , 
h{9) = p- 1 {l-9 2 ) f}5 F(z{9)). 



Table 3 

Comparison of determinations of universal ra- 
tios of amplitudes of the three-dimensional Ising 
model. 

Uo Qc Uj 



IHT 0.530(3) 
HT+LT 0.523(9) 
MC 0.560(10) 
MC 0.550(12) 
FT 0.540(11) 



0.3330(10) 1.961(7) 

0.324(6) 1.96(1) 

0.328(5) 1.95(2) 

0.331(9) 2.013(28) 



BM 
MS 
LV 



0.56(2) 
0.51(3) 
0.538(17) 



0.33(5) 
0.35(4) 



1.93(7) 
1.92(15) 



In the Ising case, corresponding to the break- 
ing of a discrete symmetry, 9q is a simple zero of 
h(9). We approximate h{&) with an odd polyno- 
mial in 9, fixing its coefficients from the small- z 
expansion of F{z). p is a free parameter; as long 
as we keep the parametric representation exact 
its value is immaterial, but it becomes significant 
once we make approximations, p can be used to 
optimize the approximation, and it can be deter- 
mined from a global stationarity condition 

We use the values of /3, S, re, r$, rio obtained 
by IHT to compute successive approximations to 
h(9); we check the stability of the values of several 
universal amplitude ratios in order to select the 
best approximation. Among the many amplitude 
ratios which can be computed from h(9), we re- 
port in Table |[/ = A+/A-,Q C = B 2 {f+f/C+, 
= f + /f~; many more ratios can be found in 
Ref. @. HT+LT is a combination of HT and low- 
temperature expansion |l^,|l^l ; the other theoreti- 
cal determinations are the same discussed for the 
critical exponents, and are taken from Refs. Q 
(IHT), [|u|n| (MC), and |§P[ (FT). For exper- 
imental data, see Refs. [§ and The agreement 
between the different determinations is again sat- 
isfactory. 

In the XY case, corresponding to the breaking 
of a continuous symmetry, 9q is a double zero of 
h(9). We therefore set 

h{9) = 9(1 - 9 2 /9 2 ) 2 {l + c 2 9 2 + c 4 # 4 ...). 

We fix 9q , C2, ... from the small- z expansion of 
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Table 4 

Comparison of determinations of the universal ra- 
tio Uq = A + /A~ of the three-dimensional XY 
model. 

A+/A- 



REFERENCES 



4 He 


1.0442 


IHT+ 4 He 


1.055(3) 


IHT* 


1.062(4) 


FT 


1.056(4) 


e-exp 


1.029(13) 



F(z), and p from the requirement h(9) w (9q — 9) 2 
for 9 -> 9 . 

Only the ratio A + /A~ is measured experimen- 
tally to high precision M. We report a selection 
of theoretical determinations: IHT+ 4 He is our 
IHT computation, using as input for a the ex- 
perimental value a = -0.01285(38) IHT* 
is a complete IHT computation, without exper- 
imental input M; FT is a g expansion in fixed 
dimension |16|; e-exp is obtained by e expansion 
|l7j ]. The value of A + /A~ is strongly correlated 
with the value of a, and all disagreement between 
IHT* and experiment can be reconduced to the 
discrepancy in a. 

4. CONCLUSIONS 

The study of HT series of "improved" models, 
with parameters determined by MC simulations, 
allowed us to compute with high precision the 
universal quantities (critical exponents and effec- 
tive potential) characterizing the critical behavior 
of the symmetric phase. 

Suitable approximation schemes allow the re- 
construction of the critical equation of state start- 
ing from the symmetric phase; many universal 
amplitude ratios can be computed. 

For the Ising universality class, theoretical 
computations are much more precise than exper- 
iments. On the other hand, for the XY class, 
some very precise experimental results for a and 
A + /A~ have been obtained 0. There is disagree- 
ment with the most precise theoretical results ||] . 
A new-generation experiment is in preparation 
p0| ]; it would be interesting to improve further 
the theoretical computations as well. 
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